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The Fokker-Planck equation (FPE) is a familiar model in various problems of nonequilibrium statistical 
OO . physics Jit] . In this paper we consider the FPE of the form 
Q\ ■ 

d t W = d x -{D-[Wd x U + d x W]}. (1) 

Here W(x,t) is the probability density over the configuration space x, at the time t, while U(x) and D(x) are 
, the potential and the positively semi-definite (y ■ D ■ y > 0) diffusion matrix. The dot denotes convolution in 

the configuration space. The FPE (Q) is particularly important in studies of polymer solutions Let us recall 
| the two properties of the FPE (ffl), important to what will follow: (i). Conservation of the total probability: 

JW(x,t)dx = 1. (ii). Dissipation: The equilibrium distribution, W eq oc exp(— U), is the unique stationary 
; solution to the FPE (Q). The entropy, 

3: s[w} = - [w( X ,t)J w{x ^ 



o 



W eq {x) 



dx, (2) 



is a monotonically growing function due to the FPE and it arrives at the global maximum in the equilibrium. 
qq ■ These properties are most apparent when the FPE (|l|) is rewritten as follows: 

0\ ; 

"w) 1 d t W(x,t)=M w -^f^ ( 3 ) 

o ; 

where Mw = —d x ■ [W(x,t)D(x) ■ d x ] is a positive semi-definite symmetric operator with kernel 1. The form 
(||) (the dissipative vector field is a metric transform of the entropy gradient) is an example of the dissipative 
""O part of a structure termed GENERIC in a recent series of papers || . 

Usually one is interested in dynamics of moments of the distribution function W rather than in the dynamics 
^ , of the W itself. Except for simplest potentials U and diffusion matrices D, the moment equations, as they 
k> ' follow from the FPE (^), are not closed. Therefore, closure procedures are required. 

In this paper we address the problem of closure for the FPE (|j) in a general setting. First, we review the 
maximum entropy principle (MEP) as a source of suitable initial approximations for the closures. We also 
discuss a version of the MEP, valid for a near-equilibrium dynamics, and which results in explicit formulae for 
arbitrary U and D. 

The MEP closures are almost never invariants of the true moment dynamics, and corrections to the MEP 
closures is the central issue of this paper. For this purpose, we apply the method of invariant manifold 
which is carried out (subject to certain approximations explained below) to explicit recurrence formulae for 
one-moment near-equilibrium closures for arbitrary U and D. These formulae give a method for computing 
the lowest eigenvalue of the problem, and which dominates the near-equilibrium FPE dynamics. 

MEP closures Let M = {Mo, Mi, . . . , M^} be linearly independent moments of interest, Mi[W] = 
J rrii (x) W{x)dx, and where mo = 1. We assume existence a function W*(M,x) which extremizes the entropy 
S (0) under the constraints of fixed M. This MEP distribution function may be written 



W* = W eq exp 



hirrii{x) — 1 



.4=0 



where A = {Ao, Ai, . . . , A^} are Lagrange multipliers. Closed equations for moments M are derived in two 
steps. First, the MEP distribution is substituted into the FPE (0) or (|J) to give a formal expression: dtW* = 
Mw* (SS/SW) | w _ w , • Second, applying a projector II*, 



IT. = ^2(dW*/dMi) [ mi(x) • dx, 
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on both sides of this formal expression, we derive closed equations for M in the MEP approximation. Further 
processing requires an explicit solution to the constraints, J W*(A, x)mi(x)dx = Mi, to get the dependence of 
Lagrange multipliers A on the moments M. Though typically the functions A(M) are not known explicitly, 
one general remark about the moment equations is readily available. Specifically, the moment equations in the 
MEP approximation have the form: 

.7 dS*(M) 
j=0 3 

where S*(M) = S[W*(M)] is the macroscopic entropy, and where M*- is an M-dependent (fc + 1) X (k + 1) 
matrix: 

m = I ' W*{M,x)[d x mi{x)] ■ D{x) ■ [d x mj{x)]dx. 



The matrix M*j is symmetric, positive semi-definite, and its kernel is the vector 5oi. Thus, the MEP closure 
reproduces the GENERIC structure on the macroscopic level, the vector field of macroscopic equations (|]) is a 
metric transform of the gradient of the macroscopic entropy. 

Triangle MEP closures || . The following version of the MEP makes it possible to derive more explicit results 
in a general setting: In many cases, one can split the set of moments M in two parts, Mj — {Mq, Mi, . . . , Mi} 
and Ma = {Af;+i, . . . , Mf.}, in such a way that the MEP distribution can be constructed explicitly for Mj 
as W?(Mi,x), The full MEP problem for M = {M/,M//} in the "shifted" formulation reads: extremize the 
functional S[Wj + AW] with respect to AW, subject to the constraints Mi[Wf + AW] = Mj and M H [Wf + 
AW] — Ma. Let us denote as AMn deviations of the moments Mu from their values in the state Wf . For 
small deviations, the entropy is well approximated with a quadratic functional AS 1 [AW] which is an expansion 
of the functional (§) in the state Wj up to the terms of the order AW 2 . With Mi[Wj] = Mi, we come 
to the following problem: extremize the functional A5[AW], subject to the constraints M/[AW] = 0, and 
Mj/[AW] = AM/j. The solution to the latter problem is always explicitly found from a (k + 1) x (k + 1) system 
of linear algebraic equations for Lagrange multipliers. 

In the remainder of this paper we deal with one-moment near-equilibrium closures: Mi = M , (i. e. Wf = 
W eq ), and the set Mu contains a single moment M — J mWdx, m(x) ^ 1. We will specify notations for the 
near-equilibrium FPE , writing the distribution function as W = W eq (l + \&), where the function 4" satisfies an 
equation: 

= w- q \m, (5) 

where J = d x ■ [W eq D ■ d x ]. The triangle one-moment MEP function reads: 



= W eq 



1 + AMto (0) 



(6) 



where AM — M — (m) , and 

mS^ — [(mm) — (m) 2 ] _1 [m— (to)]. (7) 

Brackets (...} = j W eq ... da; denote equilibrium averaging. The superscript (0) indicates that the triangle MEP 
function (1q) will be considered as an initial approximation to a procedure which we address below. Projector 
for the approximation (^|) has the form 



rt Q) * = W eq / 2 Trm l^.dx. (8) 



' eq (m(°)m(°)) 

substituting the function (JsJ) into the FPE (|J), and applying the projector (||) on both the sides of the resulting 
formal expression, we derive an equation for M: M = — XqAM, where 1/Ao is the inverse effective time of 
relaxation of the moment M to its equilibrium value, in the MEP approximation (|^): 

A = (fflOmt 01 )- 1 ^ • D ■ d x m^). (9) 

Invariant closures. Both the MEP and the triangle MEP closures are almost never invariants of the FPE 
dynamics. That is, the moments M of solutions to the FPE (|l]) vary in time differently from the solutions to the 
closed moment equations like (Q), and these variations are generally significant even for the near-equilibrium 
dynamics. Therefore, we ask for corrections to the MEP closures to finish with the invariant closures 0. 

First, the invariant one-moment closure is given by an unknown distribution function W^ 00 ^ = W eq [l + 
AMm' 00 ' (x)] which satisfies an equation 
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[l-n< <x0 ]J"* (oc) =0- (10) 

Here n(°°) is a projector, associated with an unknown function m(°°\ and which is also yet unknown. Eq. ( |10| ) 
is a formal expression of the invariance principle for a one-moment near-equilibrium closure: considering W^> 
as a manifold in the space of distribution functions, parameterized with the values of the moment M , we require 
that the microscopic vector field JrrS 00 ^ be equal to its projection, n.(°°> Jm^°°\ onto the tangent space of the 
manifold W^°°\ 

Now we turn our attention to solving the invariance equation (|l0|) iteratively, beginning with the triangle 
one-moment MEP approximation (||). We apply the following iteration process to the Eq. (|l0|): 

[i-n«]jm( fe+1 ) = 0, (11) 

where k = 0, 1, . . ., and where m^ k+1 ' = rrS k ^ + ^- k+1 \ and the correction satisfies the condition {^- k+1 ^ [ rrS k ^) = 0. 
Projector is updated after each iteration, and it has the form 



n(fc+1) « - W « {m W m W) J m(k+1){x) ' dX - (12) 

Applying n( fc+1 ) to the formal expression, W / eg r7i/ fe+1 ) M = AM[1 — n^ fc+1 ^] JrrS k+1 \ we derive the macroscopic 
equation, M ~ — Afc+iAM, where A^+i is the (k + l)th update of the inverse effective time (||): 

_ (d x m^ ■ p-d.M k+1) ) n „. 

Specializing to the one-moment near-equilibrium closures, and following a general argument Qj, solutions to 
the invariance equation ( |l0| ) are eigenfunctions of the operator J, while the formal limit of the iteration process 
( |TT| ) is the eigenfunction which corresponds to the eigenvalue with the minimal nonzero absolute value. 

Diagonal approximation. To obtain more explicit results in the iteration process (Ul|), we introduce an 
approximate solution on each iteration. The correction ^S k+1 ^ satisfies the condition (rn^^/l ) = 0, and can 
be decomposed as follows: ^ k+r > = a k e^ + e (k } t . Here = W^\\ - Tl^Jm^ = X k m^ + R.W is the 
variance of the fcth approximation, where 

= d x ■ [D ■ d x m {k) ] - d x U ■ D ■ d x m {k) . (14) 

The function e^} t is orthogonal to both e( fc ) and m^: (e^e^) — 0, and (mi k ^ e^ r \) — 0. Our diagonal 
approximation (DA) consists in disregarding the part e^ t . Specifically, we consider the following ansatz at the 
kth iteration: 

m (fc+i) =m W+a fe e (fe) . (15) 

Substituting the ansatz (15) into the Eq. jTl|), and integrating the latter expression with the function e^ k \ we 
evaluate the coefficient a k : 

ak = A3-1Ui*+B fc ' (16) 
where parameters A k and B k represent the following equilibrium averages: 

A k = (mC'mW)- 1 ^^') (17) 
B k = {m ik) m {k) }- 1 (d x R (k) ■ D ■ d x R {k) ). 

Finally, putting together Eqs. @, (Q), ©, ©, and (0), we arrive at the following DA recurrency 
solution, and which is our main result: 

m (*+i) _ m W + ak [X km ^ + , (18a) 
_ A fc - {A k - X 2 k )a k 
Xk+1 l + (A k - Xl)al ■ (18b) 

To test the convergency of the DA process ( |l8| ) we have considered two potentials U in the FPE (Q) with a 
constant diffusion matrix D. The first test was with the square potential U = x 2 /2, in the three-dimensional 
configuration space, since for this potential the detail structure of the spectrum is well known. We have 
considered two examples of initial one-moment MEP closures with mA°) = X\ + 100(x 2 — 3) (example 1), and 



(k+i) 
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to' ) — x\ + lOOrr 6 ^ (example 2), in the Eq. (Q). The result of performance of the DA for (18b) is presented 
in the Table Q, together with the error Sf. which was estimated as the norm of the variance at each iteration: 
$k = {e^e^)/{m^m^}. In both examples, we see a good monotonic convergency to the minimal eigenvalue 
Aoo = 1) corresponding to the eigenfunction x\. This convergency is even striking in the example 1, where the 
initial choice was very close to a different eigenfunction x 2 — 3, and which can be seen in the non-monotonic 
behavior of the variance. Thus, we have an example to trust the DA as converging to the stationary point of 
the original iteration procedure (pi]). 

For the second test, we have taken a one-dimensional potential U = — 501n(l — x 2 ), the configuration 
space is the segment \x\ < 1. Potentials of this type (so-called FENE potential) are used in applications of 
the FPE to models of polymer solutions Q . Results are given in the Table |l| for the two initial functions, 
TO (o) _ x 2 _|_ iQ x i _ (^2 _|_ i Ox 4 } (example 3), and = x 2 + 10x 8 — (x 2 + 10a; 8 ) (example 4). Both the examples 
demonstrate a stabilization of the Xk at the same value after some ten iterations. 

In conclusion, we have developed the principle of invariance to obtain moment closures for the Fokker-Planck 
equation (Q), and have derived explicit results for the one-moment near-equilibrium closures, particularly 
important to get information about the spectrum of the FP operator. 
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TABLE I. Iterations X k and the error 8 k for U = x 2 /2. 
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12 


16 


20 


Ex. 1 


A 


1.99998 


1.99993 


1.99575 


1.47795 


1.00356 


1.00001 


1.00000 




6 


0.16 ■ 10~ 4 


0.66 ■ 10~ 4 


0.42 ■ 


0.24 


0.35 ■ 10"* 


0.13 ■ 10~ 4 


0.54- 10~' 
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Ex. 2 


X 


3.399 


2.437 


1.586 


1.088 


1.010 


1.001 


1.0002 




8 


1.99 


1.42 


0.83 


0.16 


0.29 ■ 10~ 4 


0.27 ■ 10~^ 


0.57 ■ W- 3 



TABLE II. Iterations A fe for U = -501n(l - x 2 ). 
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Ex. 3 


A 


213.1774 


212.1864 


211.9148 


211.8619 


211.8499 


211.8453 


211.8433 


211.8422 


211.8417 


Ex. 4 


A 


216.5856 


213.1350 


212.2123 


211.9984 


211.9295 


211.8989 


211.8838 


211.8757 


211.8713 
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